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Abstract A gravitational wave stochastic background of astrophysical origin may have 
resulted from the superposition of a large number of unresolved sources since the be- 
ginning of stellar activity. Its detection would put very strong constrains on the physical 
properties of compact objects, the initial mass function or the star formation history. On 
the other hand, it could be a 'noise' that would mask the stochastic background of cosmo- 
logical origin. We review the main astrophysical processes able to produce a stochastic 
background and discuss how it may differ from the primordial contribution by its statisti- 
cal properties. Current detection methods are also presented. 
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1 INTRODUCTION 

Gravitational wave (GW) astronomy will enable a new window to the Universe to be opened: not only 
one expects to discover a set of new exotic sources, but also to travel back in time, toward the very early 
stages of the evolution of the Universe. According to various cosmological scenarios, we are bathed in 
a stochastic background of gravitational waves, memory of the first instant of the Universe, up to the 
limits of the Plank era and the Big Bang (Grishchuk et al. 2001). Proposed theoretical models include 
the amplification of vacuum fluctuations during inflation (Grishchuk 1974, 1993; Starobinskii 1979), pre 
Big Bang models (Gasperini & Veneziano 1993, 2003; Buonanno et al. 1997), cosmic strings (Vilenkin 
& Shellard 2000) or phase transitions Caprini (2010) (see Maggiore 2000 or for a general review and 
references therein). In addition to this cosmological background (CGB), an astrophysical contribution 
(AGB) may have resulted from the superposition of a large number of unresolved sources since the 
beginning of the stellar activity, which can be either short live burst sources, such as core collapses to 
neutron stars (Blair & Ju 1996; Coward et al. 2001, 2002; Howell et al. 2004; Buonanno et al. 2005; 
Marassi et al. 2009; Zhu, Howell & Blair 2010) or black holes (Ferrari et al. 1999a; de Araujo et al. 
2000, 2002a,b, 2004), oscillation modes (Owen et al. 1998; Ferrari et al. 1999b; Marassi et al. 2009; 
Zhu, Howell & Blair 2010), final stage of compact binary mergers (Regimbau & de Freitas Pacheco 
2006b; Regimbau 2007), or periodic long live sources, typically pulsars (Regimbau & de Freitas Pacheco 
2001a, 2006a), the early inspiral phase of compact binaries (Ignatiev et al. 2001; Schneider et al. 2001; 
Farmer & Phinney 2002; Cooray 2004) or captures by supermassive black holes (Barack & Cutler 2004; 
Schnittman et al. 2006), whose frequency is expected to evolve very slowly compared to the observation 
time. The nature of the AGB may differ from its cosmological counterpart, expected to be stationary, 
unpolarized, gaussian and isotropic, by analogy with the cosmic microwave background (CMB). On 
the one hand the distribution of galaxies up to 100 Mpc is not isotropic but strongly concentrated in 
the direction of the VIRGO cluster and the Great attractor; on the other hand, depending whether the 
time interval between events is short compared to the duration of a single event, the integrated signal 
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may result in a continuous, a popcorn noise or a shot noise background (Coward & Regimbau 2006). 
The optimal strategy to search for a Gaussian (or continuous) stochastic background is to cross corre- 
late measurements of multiple detectors, which can be either resonant antennas such as the cryogenic 
bars AURIGA, NAUTILUS, EXPLORER, ALLEGRO or NIOBE (Amaldi et al. 1990; Cerdonio et al. 
1997; Pallottino 1997; Mauceli et al. 1996; Blair et al. 1995), laser interferometers such as LIGO, 
VIRGO, GEO600, TAMA300/LCGT and the third generation Einstein Telescope (Abramovici et al. 
1992; Bradaschia et al. 1990; Hough 1992; Kuroda et al. 2006) on earth or LISA in space (Bender & the 
LISA Study Team 1998), or natural detectors such as millisecond pulsars of the Parkes Pulsar Timing 
Array (PPTA) (Jenet et al. 2005; Manchester 2006). Space and terrestrial detectors will be complemen- 
tary in the 10 -5 — 10 4 Hz band, while the PPTA is expected to detect GWs at nHz frequencies. Over the 
last decade, the first generation of terrestrial detectors have been built, commissioned and are running 
in scientific mode at (or close to) their design sensitivities, providing the opportunity to do joint data 
analysis. 

This chapter gives an overview of the main features of the AGBs, from modeling to detection, 
and discuss how their statistical properties may differ from those of the cosmological background. In 
Section 1, we review the spectral and statistical properties of astrophysical backgrounds, in Section 2 we 
introduce the actual detection method to search for a stochastic background in a network of detectors, 
in Section 3 we review the most popular predictions of the cosmological background, in Section 4 
we describe models of AGB and in Section 5 we discuss the current observational results and future 
prospects. 



2 THE CHARACTERISTICS OF THE SPECTRUM 

The spectrum of the gravitational stochastic background is usually characterized by the dimensionless 
parameter (Allen & Romano 1999): 

o f..\ _ 1 d Ps™ (1) 



p c d In v Q 

where p gw is the gravitational energy density, v a the frequency in the observer frame and p c — 

the critical energy density needed to make the Universe flat today. Experimentalists may prefer to work 

with the spectral energy density, 

2>H 2 1 

which is directly comparable to the detector sensitivity. 

For a stochastic background of astrophysical origin, the energy density parameter is given by 
(Ferrari et al. 1999a): 

where the integrate flux received on Earth (in erg/Hz _1 //cm -2 ), at the observed frequency v is defined 

as: 

Fvo = / p(8)U o (d,z^o) dR °f z ,Z) d0dz, (4) 

where p(6) is the probability distribution of the source parameters 9. The first factor in the integral is 
the fluence of a source (in erg/Hz^Vcm -2 ) located at redshift z: 

1 rl 7? 

f^ z ^ = ^ r{z) 2 d ; w (^ (i+*)), (5) 

where r (z) is the proper distance, which depends on the adopted cosmology, dE ^J (6, v) the gravitational 
spectral energy emitted and v a (1 + z) the frequency in the source frame. The second factor is the number 
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of sources in the interval 9 — 9 + d9, per unit of time in the observer frame and per redshift interval, is 
given by: 

dR°(9,z) dV 

___ = „( M )_( z)> (6) 

where p°(9, z) the event rate in Mpc~ 3 yr _1 and ^r{z) the comoving volume element. 
Combining the expressions above, one obtains for the density parameter: 



<W*>) = j dO P (9) to^-^QMl + z)). (7) 

Replacing the constants by their usual values we get: 



r! g wK) = 5.7x 10- 5b (-^ Yv J d9p(9)j dz ^ ^ ^ {6, z/ (l + z)), (8) 

where p° is given for ho = 0.7. The limits of the integral over z depend on both the emission frequency 
range in the source frame, and the redshift interval, where the source can be located: 

if v < 



z (9 v ) = i (i+w) (a) 
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and 

if v > 



Consequently, the shape of the spectrum of any astrophysical background is characterized by a cutoff 
at the maximal emission frequency and a maximum at a frequency which depends on the shape of both 
the redshift distribution and the spectral energy density. 

For most of the models presented in Section 5, the event rate per unit of redshift can be derived 
directly from the cosmic star formation rate. In the simple case when the gravitational emission occurs 
shortly after the birth of the progenitor, it is given by: 

p°{0,z)=\(O,z)£&, (11) 
1 + z 

where A is the mass fraction converted into the progenitors in Mq 1 , which depends on the initial 
mass function, ^ the element of comoving volume and p* (z) the cosmic star formation rate (SFR) 
in M Mpc~ 3 yr -1 . The (1 + z) factor in the denominator corrects for the time dilatation due to the 
cosmic expansion. 

Observations of star forming galaxies with large telescopes such as the Keck or the Hubble Space 
Telescope have extended our view of the Universe up to redshifts z ~ 5 — 6, by tracing the evolution 
with cosmic time of the galaxy luminosity density. The main uncertainty comes from dust extinction, 
which spreads the UV luminosity into the far IR. Madau et al. (1998) derived an expression that matches 
most of the measurements in the U-V continuum and Ha, up to z ~ 4, and that includes an extinction 
correction of Ai5oo=l-2 mag. The SFR is expected to increase rapidly between z ~ — 1, peak around 
z ~ 1.7 and smoothly decrease at large redshifts. After z ~ 2, the behavior must be regarded as tentative, 
due to the large uncertainties in the estimates of the U-V luminosity from the Lyman break galaxies in 
the Hubble Deep Field. Steidel et al. (1999) proposed another scenario where the SFR remains constant 
after z ~ 2. Other studies suggested even an increase of the SFR, claiming that it may have been severely 
underestimated due to large amount of dust extinction (Blain et al. 1999). However, the hypothesis of 
a gentle decline at high redshifts seems to be favored by new measurements of the galaxy luminosity 
function in the UV (SDSS, GALEX, COMBO 17) and FIR wavelengths (Spitzer Space Telescope), 
which allowed to refine the previous models of star formation history, up to redshift z ~ 6, with tight 
constraints at redshifts z < 1. In a recent work, Hopkins & Beacom (2006) used the Super Kamiokande 
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limit on the electron antineutrino flux from past core-collapse supernovas to derive parametric fits of the 
form of Cole et al. (2001). Investigating the effect of the initial mass function (IMF) on the normalization 
of the SFR, they showed that top heavy IMFs are preferred to the traditional Salpeter IMF (Salpeter 
1955), and the fits are optimized for IMFs of the form: 



with a turnover below too = 0.5 M , normalized within the mass interval 0.1 — 100 M Q such as 
/ m£(m)dm = 1, and with 7 = 2.35 (A modified Salpeter). Fardal et al. (2007) used a different set of 
measurements and different dust extinction corrections and found an SFR similar to that of Hopkins & 
Beacom (2006) up to z ~ 1, but which decreases slightly at higher redshifts. Wilkins et al. (2008) used 
measurements of the stellar mass density and derived an SFR equivalent to that of Hopkins & Beacom 
(2006); Fardal et al. (2007) for redshifts smaller than z ~ 0.7, but again is lower at higher redshifts. 
Finally, Nagamine et al. (2006), derived a model from the fossil record of star formation in nearby 
galaxies. It is probably underestimated at small redshifts, and is constant at high redshifts due to the 
contribution of elliptical galaxies. Note that at present there is a discrepancy between the "instantaneous" 
SFR, measured from the emission of young stars in star forming regions, and the SFR as determined 
from extragalactic background light. This could have an important impact on the contribution to the 
confusion background for sources from z > 2. However, it shouldn't noticeably affect the results, since 
sources beyond z ~ 2 are too weak to contribute significantly to the integrated signal. Figure compares 
the four prior models described above, calculated for the flat Einstein de Sitter 737 cosmology, with 
fi m = 0.3, — 0.7 and Hubble parameter Hq = 70 km s _1 Mpc -1 (Rao et al. 2006), corresponding 
to the so-called concordant model derived from observations of distant type la supernovae (Perlmutter 
et al. 1999) and the power spectra of the cosmic microwave background fluctuations (Spergel et al. 
2003). 

Besides the spectral properties, it is important to study the nature of the background. In the case of 
short-lived signals, they may show very different statistical behavior depending on the ratio between the 
duration of the events and the time interval between successive events, the duty cycle: 



which is also the average number of events present at the detector at a given observation time. 

Continuous: the number of sources is large enough for the time interval between events to be small 
compared to the duration of a single event. The waveforms overlap to create a continuous background 
and due to the central limit theorem, such backgrounds obey the Gaussian statistic. They are completely 
determined by their spectral properties and could be detected by data analysis methods in the frequency 
domain such as the standard cross correlation statistic (Allen & Romano 1999). 

Shot noise: the number of sources is small enough for the time interval between events to be long 
compared to the duration of a single event. The waveforms are separated by long stretches of silence 
and the closest sources may be detected by data analysis techniques in the time domain (or the time 
frequency domain) such as match filtering (Arnaud et al. 1999; Pradier et al. 2001). 

Popcorn: an interesting intermediate case arises when the time interval between events is of the 
same order of the duration of a single event. These signals, which sound like crackling popcorn, are 
known as "popcorn noise". The waveforms may overlap but the statistic is not Gaussian anymore and the 
amplitude on the detector at a given time is unpredictable. Promising data analysis strategies have been 
investigated in the last few years, such as the Maximum Likelihood statistic, an extension of the cross 
correlation statistic in the time domain (Drasco & Flanagan 2003) or methods based on the Probability 
Event Horizon concept (Coward & Burman 2005), which describes the evolution of the cumulated signal 
throughout the Universe, as a function of the observation time. 




for 0.1 < to < to 
for toq < to < 100 



(12) 




(13) 
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Fig. 1 Cosmic star formation rates (in M Mpc 3 yr ) used in this paper: Hopkins & 
Beacom (2006) (continuous line), Fardal et al. (2007) (dashed line), Wilkins et al. (2008) 
(dot-dashed line), and the fossil model of Nagamine et al. (2006) (dot line). As discussed in 
the text, these rates are largely the same up to z ~ 1, but show important differences at higher 
redshift. 



3 DETECTION 

The optimal strategy to search for a gaussian (or continuous) stochastic background, which can be 
confounded with the intrinsic noise background of the instrument, is to cross correlate measurements 
of multiple detectors. In this section, we give a brief overview of the standard data analysis technique 
currently used for terrestrial interferometers. 

When the background is assumed to be isotropic, unpolarized and stationary, the cross correlation 
product is given by Allen & Romano (1999): 



Y= / sV(f)Q(f)s 2 (f)df, (14) 

where 

Q(f) oc (15 ) 

is a filter that maximizes the signal to noise ratio (S/R). In the above equation, P\{f) and P2(f) are 
the power spectral noise densities of the two detectors and T is the non-normalized overlap reduction 
function, characterizing the loss of sensitivity due to the separation and the relative orientation of the 
detectors (Fig. , The optimized S/N ratio for an integration time T is given by Allen (1997): 
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Fig. 2 time series corresponding to shot noise, popcorn and gaussian regimes 



In the literature, the sensitivity of a pair of detectors is usually given in terms of the minimum 
detectable amplitude corresponding to Q gw equal to constant (hereafter flat spectrum) (Allen & Romano 
1999): 

4^ 2 . _ n ,_ z 100 r 2 (/) 



3^o 2 



(erfc _i (2/3) - erfc _1 (2a))[ / df 



-1/2 



(17) 



The expected minimum detectable amplitudes for the main terrestrial interferometer pairs, at design 
sensitivity (Fig. , and after one year of integration, are given in Table , for a detection rate a = 90% 
and a false alarm rate (3 = 10%. 

^min is of the order of 10~ 6 — 10 -5 for the first generation of interferometers combined as 
LIGO/LIGO and LIGO/Virgo. Their advanced counterparts will permit an increase of two or even three 
orders of magnitude in sensitivity (fi m i n ~ 10~ 9 — 10 -8 ). The pair formed by the co-located and co- 
aligned LIGO Hanford detectors, for which the overlap reduction function is equal to one, is potentially 
one order of magnitude more sensitive than the Hanford/Livingston pair, provided that instrumental and 
environmental noises can be removed. In Table we show the evolution of the upper limit obtained 
with the LIGO detectors in a narrow band around 100 Hz, and corresponding to £! gw equal to constant 
at all frequencies. In Table , the latest published LIGO upper limit is compared to observational limits 
already achieved with resonant bar experiments at about 900 Hz and pulsar timing at nHz frequencies . 

An extension of the cross-correlation method to non isotropic contributions has been investigated 
by Allen & Ottewill (1997); Cornish (2001); Ballmer (2006) and Mitra et al. (2008). The basic idea is to 
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Fig. 3 overlap reduction function for the most promising detector pairs. L stands for LIGO 
Livingston and H for LIGO Hanford, V for Virgo, G for GEO600 and ET for the planned 
Einstein Telescope in the triangular configuration. 
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Fig. 4 designed sensitivities of the main first generation interferometers, compared to the 
planned sensitivities of advanced detectors (LIGO or Virgo) and third generation detectors 
(Einstein Telescope). 
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Table 1 Expected il m - m for the main detector pairs, corresponding to a fiat background spec- 
trum, one year of integration over the full frequency band, a detection rate a = 90% and a 
false alarm rate (3 = 10%. LHO and LLO stand for LIGO Hanford Observatory and LIGO 
Livingston Observatory, ET stand for ET for the planned Einstein Telescope in the triangular 
configuration (Flanagan 1993). 
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Table 2 evolution of the LIGO 90% Bayesian upper limit on a frequency independent fi, 



Run 


frequency band (Hz) 


upper limit 


reference 


SI 


40-314 


23 


Abbott et al. (2004) 


S3 


69 - 156 


8.4 x 10" 4 


Abbott et al. (2005) 
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51 - 150 


6.5 x 10~ 5 


Abbott et al. (2007) 
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40 - 170 


5.9 x 10" 6 


Abbott et al. (2009) 



Table 3 best published direct upper limits on a frequency independent fl gw derived from 
correlation, for different type of experiments. 



Type of detectors 


Experiment 


frequency (Hz) 


upper limit 


reference 


Room Temp, resonant bars 


Glasgow 


985 


6125 


Hough etal. (1975) 


Cryogenic resonant bar 


Explorer+Nautilus 


907 


120 


Astone etal. (1999) 


Pulsar timing 


Parkes 


1(J -8 


4 x 10~ 8 


Jenet et al. (2005) 



use multiple detector pairs to create maps of anisotropy of the GW background, similar to a radiometer 
ofGWs. 



4 RELIC STOCHASTIC BACKGROUND 

Mechanisms able to generate stochastic backgrounds of GWs in the very early stages of the Universe 
have been investigating intensively in the past decades. Their detection would have a profound impact 
on our understanding of near Big Bang cosmology and high energy physic, providing a unique way 
to explore the Universe a fraction of second after the Big Bang, after gravitons decouple from the 
primordial plasma. It is not the purpose of this article to develop in details all the different models of 
cosmological stochastic background present in the literature, as our main interest is the astrophysical 
background, but for comparison purpose, we give in this section a rapid overview of some popular 
predictions that could be masked by the astrophysical background. We refer interested readers to very 
nice review papers by Allen (1997); Maggiore (2000) and Buonanno (2003) . 

In this section, unless it is mentioned otherwise, the Hubble parameter is assumed to be Hq = 
70kms _1 Mpc~l. 

4.1 Amplification of Vacuum Fluctuations 

Amplifications of vacuum fluctuations at the transitions between the de Sitter, radiation dominated (RD) 
and matter dominated (MD) eras, first discussed by Grishchuk (1974, 1993) and Starobinskii (1979), are 
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Fig. 5 Theoretical predictions of the cosmological stochastic background and observational 
bounds. The cosmic strings plot corresponds to p = 0.1, £ = 7 x 10~ 5 , G[i = 10~ 8 . This 
figure was kindly provided by Vuk Mandic. 



expected to produce a GW background which spectrum and amplitude depend strongly on the fluctuation 
power spectrum developed during the early inflationary period. The standard de Sitter inflation predicts 
a spectrum that decreases as l/f 2 in the range 3 x 10~ 18 — 10 -16 Hz and then remains constant in a very 
large band up to MHz frequencies. In the low frequency region, modes amplified during inflation at both 
the de Sitter-RD and RD-MD transitions contribute. The turnover between the two phases corresponds 
to the limit at which only the modes amplified during the de Sitter-RD transition can be observed. The 
COBE experiment, which has the same l/f 2 behaviour in the low frequency region, provides an upper 
bound of fJ gw - 9 x 10" 14 (Maggiore 2000; Buonanno 2003) for the flat region. Actually, a GW 

background larger than f2 gw ~ 7x 10~ 10 (^j — ) 2 at frequencies between 3 x 10~ 18 — 10~ 16 Hz would 
have produced stochastic frequency redshifts through Sachs-Wolfe effects, which would have resulted 
into temperature fluctuations larger than those measured for the cosmic microwave background. In the 
more realistic scenario of "slow roll down" inflation, the inflaton field rolls toward the minimum of 
its potential, producing an acceleration of the expansion. The Hubble parameter is not constant like in 
the standard scenario, but decreases monotonically during the period of inflation. GWs are produced by 
fluctuations that go out the Hubble radius during inflation, and re-enter at the radiation era. The resulting 
spectrum is not flat as for the de Sitter inflation but rather has a f lT dependency, where tit < and 
\tit\ << 1. The spectral index can be expressed in terms of the scalar and tensorial contributions to 
the quadrupole cosmic microwave background (CMB) anisotropy as nr = —T/7S (Maggiore 2000). 
The most optimistic predictions for a detection with LISA at / ~ 10~ 4 , corresponding to rij< = 0.175, 
give an amplitude £l gw ~ 10~ 15 , but tit could be much smaller, of the order of 10 -3 (Maggiore 2000). 
In a recent paper, Tong & Zhang (2009) studied the effect of a running spectral index at on the GW 
spectrum and found that a t > could enhance the signal significantly, especially at high frequencies. 
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A more interesting case arises from pre-big-bang scenarios in string cosmology (Gasperini & 
Veneziano 1993, 2003). According to these models, the standard RD and MD eras were preceded by 
phases in which the Universe was first large and shrinking (inflaton phase) and then characterized by a 
high curvature (stingy phase). The GW spectrum produced at the transition between the stingy phase and 
the RD era is described as 51 gw ~ f 3 for f < f s and f2 gw ~ j3-2^i f or / s < / < / x (Buonanno et al. 
1997; Mandic & Buonanno 2006). The turnover frequency is essentially unconstrained, /i < 1.5 reflects 
the evolution of the Universe during the 'stingy' phase and the cutoff frequency /i, which depends on 
string related parameters, has a typical value of 4.3 x 10 10 Hz. An upper limit on f2 gw is imposed by 
the Big Bang Nucleosynthesis (BBN) bound down to 10~ 10 Hz, corresponding to the horizon size at the 
time of BBN. Actually, if the total energy amount carried by GWs, J fl gw d(lnf), at the time of nucle- 
osynthesis was larger than 1.1 x 10 _5 (A r ly — 3), where N v is the effective relative number of species 
at the time of BBN, it would have resulted into a particle production rate too large compared to the 
expansion of the Universe, to account for the primordial abundances of the light elements 2 H, 2 He, 4 He 
and 7 Li. Measurements of the light element abundances combined with the WMAP data gives N v < 4.4 
(Cyburt et al. 2005), which translates to Sl gw < 1.5 x 10 -5 . Recent measurements of CMB anisotropy 
spectrum, galaxy power spectrum and of the Lyman-a forest give a bound of similar amplitude which 
extends down to 10~ 15 Hz, corresponding to the horizon size at the time of CMB decoupling (Smith 
etal. 2006). 

4.2 Cosmic Strings 

Cosmic strings, formed as linear topological defects during symmetry breaking phase transitions or in 
string theory inspired inflation scenarios, may emit GWs by oscillating relativistically and shrinking 
in size (Buonanno 2003). CMB observations are not consistent with the most promising scenario of 
very large mass-per-unit-length strings, acting as initial seeds for the formation of large-scale structures 
at the GUT scale symmetry break, but strings of lower energy scale may still contribute to the CGB. 
Also, in models with non vanishing cosmological constant, it can still be a viable option (Avelino et al. 
1998; Battye et al. 1998). The spectrum is expected to peak around the frequency / ~ 10~ 12 and 
become almost fiat in a large frequency band from 10~ 8 to 10 10 Hz, where the amplitude can reach 
f! gw ~ 10~ 9 — 10~ 8 , according to numerical simulations of a cosmic strings network with G/i < 
10~ 6 (Buonanno 2003). At the present time, the most stringent constrain is given by pulsar timing 
observations. When passing between Earth and pulsars, GWs may cause fluctuations in the time arrival 
of the pulses. Observations of PSR B 1805+09 (Kaspi et al. 1994; Lommen et al. 2003), a very stable 
narrow profile millisecond pulsar, with the Arecibo and Green Bank radiotelescopes for over 17 years, 
give an upper limit of fl gw ~ 1.2x 10~ 9 at the frequency / = = 1.86 x 10 -9 nHz, and combining 
the timings of seven pulsars citepjen05 placed a lower bound of 57 gw ~ 4 x 10~ 8 . The Parkes Pulsar 
Timing Array project (Manchester 2006) which is expected to reach O gw ~2x 10~ 10 , by monitoring 
twenty pulsars for five years, may be our best hope to detect cosmic strings in the near future. 

In a recent work, Damour & Vilenkin (2000, 2001, 2005) and Siemens et al. (2007) considered 
the stochastic background created by cusps of oscillating cosmic superstring loops at the end of Brane 
inflation. The amplitude and the shape of the GW spectrum is determined by three parameters (Damour 
& Vilenkin 2005): the string tension /i, the reconnection probability p, typically in the range 10~ 3 — 1 and 
e the typical size of the closed loops produced in the string network. In particular, /i and e determine the 
lowest frequency at which a string loop can emit GWs. The GW spectrum is characterized by a decrease 
at lowest frequencies, followed by a flat region. Assuming p = 5x 10~ 3 , G\i = IQ- 7 and e = 10~ 7 , 
(Siemens et al. 2007) obtained a spectrum that avoids the low frequency bound due to CMB or pulsar 
timing measurements but still remain in the sensitivity band of space or ground based detectors. Let's 
mention that according to Damour & Vilenkin (2000, 2001) and Damour & Vilenkin (2005) occasional 
strong beams of GWs could be produced at cusps, forming a popcorn like noise on top of the gaussian 
contribution. 
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4.3 Phase Transitions 

At the early stages of its evolution, the Universe may have undergone several episodes of phase tran- 
sition, in which the symmetry of particle-physics fundamental interactions spontaneously broke. This 
may occur for instance at the QCD (150 Mev) and electroweak scales (100 GeV) or even earlier, at the 
grand unified scale (see Maggiore 2000 and references therein). The Standard Model predicts rather a 
smooth crossover, but in its supersymmetric extensions, the transition from a metastable phase (the false 
vacuum) to the state of broken symmetry (the true vacuum) can be first order, and large amount of GWs 
could be produced when bubbles of the new phase are nucleated, grow and as they become more numer- 
ous, collide at very high velocities. The GW spectrum reach a maximum of f2 gw ~ 10~ 6 (^) 2 (^) 1 / 3 

at /max ~ 3x ( if")(i^o ) 1 ^ 6 -^*' w here T = Toe - ' 3 * is the nucleation rate of bubbles, T* is the tempera- 
ture in GeV of the phase transition, H* the relevant Hubble parameter and the number of relativistic 
degrees of freedom. In particular, a phase transition at the electroweak scale could give a detectable 
signal of Sl gw ~ 10~ 9 at the mHz frequency, where LISA is the most sensitive. Besides the collision of 
the broken phase bubbles, other processes are expected to produce gravitational waves, such as the mag- 
netohydrodynamical turbulence in the plasma stirred by the bubble collisions, and the magnetic fields 
amplified by the magnetohydrodynamical turbulence (Caprini 2010). 

5 SOURCES OF ASTROPHYSICAL BACKGROUNDS 

Many examples of astrophysical backgrounds can be found in the literature. However, a direct compari- 
son between the different models is made difficult by the fact that they often use different cosmologies, 
SFRs, IMFs, or mass range for neutron star (NS) or black hole (BH) progenitors. In this section we 
review some of the most promising predictions, since it is impossible to cover in one chapter all the 
literature on the subject. 

5.1 Binary Neutron Stars 

Double neutron star coalescences, which may radiate about 10 53 erg in the last seconds of their inspiral 
trajectory, up to 1.4—1.6 kHz, may be the most important contribution in the frequency range of ground 
based detectors (Regimbau & de Freitas Pacheco 2006b; Regimbau & Mandic 2008). 

The coalescence rate per comoving volume p° in Equation 8, results from the convolution of the for- 
mation rate of the progenitors with the probability distribution P of the delay td between the formation 
of the progenitors and the coalescence: 



where z is the redshift at the time of the coalescence and Zf is the redshift at the time of formation of 
the binary. Population synthesis (Piran 1992; Tutukov & Yungelson 1994; Lipunov et al. 1995; Ando 
2004; de Freitas Pacheco et al. 2006; Belczynski et al. 2006; O'Shaughnessy et al. 2008) suggest that 
the delay time is well described by a probability distribution of the form: 



This broad model accounts for the wide range of merger times observed in binary pulsars and is also 
consistent with short gamma ray burst observations in both late and early type galaxies (Berger et al. 
2007). Belczynski & Kalogera (2001) and Belczynski et al. (2006) have identified a new efficient for- 
mation channel which produces a significant fraction of tight binaries with merger times in the range 
r m ~ 0.001 — 0.1 Myr, which gives a minimal delay time t ~ 20Myr, corresponding roughly to the 
time it takes for massive binaries to evolve into two neutron stars. 

The local cosmological rate at z = 0, p Q in Myr -1 Mpc~ 3 , is usually extrapolated by taking the 
product of the rate in the Milky Way (r mw in yr" 1 ) and the density of Milky- Way equivalent galaxies, 




(18) 



Pd(td) oc — with t d > t . 



(19) 



12 



T. Regimbau 



Table 4 Taken from Table 4 of Postnov & Yungelson (2006), most current estimates of the 
Galactic merger rates of NS-NSs and NS-BHs, derived from statistical studies (first row), and 
from population synthesis. The high rate obtained by Tutunov and Yungelson (1993) is due to 
the assumption that neutron stars or black holes are bom with no kick velocity, leading to an 
overestimate of the number of systems that survive the two supernovae. The low rate obtained 
by Voss and Tauris (2003) is due to the use of a different value of the parameter A, which 
measures the binding energy of the common envelop. 



statistics 

Kalogera et al. (2004) 


NS-NS 
83 (17-292) 




population synthesis 


NS-NS 


NS-BH 


Tutunov and Yungelson (1993) 


300 


20 


Lipunovetal. (1997) 


30 


2 


Potergies Zwart and Yungelson (1998) 


20 


2 


Nelemans et al. (2001) 


20 


4 


Voss and Tauris (2003) 


2 


0.6 


O'Shaughnessy et al. (2005) 


7 


1 


de Freitas Pacheco et al. (2006) 


17 




Belczinsky et al. (2007) 


10-15 


0.1 


O'Shaughnessy et al. (2008) 


30 


3 



given from measurements of the blue stellar luminosity around n mw ~ (1 — 2) x 10~ 2 Mpc~ 3 (Phinney 
1991; Kalogera et al. 2001; Kopparapu et al. 2008). The most current estimates of the NS-NS galactic 
coalescence rate are given in the range 1 — 817 Myr -1 (95% confidence intervals) for statistical studies 
which extrapolate the rates from observed galactic NS-NS (Kalogera et al. 2004), preferably between 
17 — 292 (95% confidence intervals) with a peak probability around 83 Myr -1 , and in the range 1 — 300, 
more likely around 10 — 30, for population synthesis models, which combine theoretical and observa- 
tional constraints (Table ). In the quadrupolar approximation, the GW energy spectrum emitted by a 
binary system, which inspirals in a circular orbit is given up to the last stable v max orbit by: 



at? IA (Gtt) 2/3 wim 2 _ 1/3 

dE gw dv = - —-i> i. (20) 

3 (mi + m,2) /A 

2/3 



Assuming m\ = = 1.4 M for the star masses, the energy density increases as v before it 
reaches a maximum of il gw ~ 3.5 x 10 _9 po at around 500 Hz, where p is the local rate in My -1 
Mpc~ 3 (about 0.01 times the galactic rate). This means that ET should be able to detect the background 
from binaries even for the most pessimistic predictions of the coalescence rate, down to po ~ 0.02 
(roughly equivalent to a galactic rate of 2 My -1 ), for a signal-to-noise ratio of 3, after one year of 
observation. 



5.2 Rotating Neutron Stars: Tri-axial Emission 

Rotating neutron stars with a triaxial shape may have a time varying quadrupole moment and hence ra- 
diate GWs at twice the rotational frequency. The total spectral gravitational energy emitted by a neutron 
star born with a rotational period Pq, and which decelerates through magnetic dipole torques and GW 
emission, is given by: 

^ = Ku 3 (l + -^v 2 )- 1 with v e [0 - 2/P ], (21) 

where 

5c 5 R e B 2 v ; 

R is the radius of the star, e — (I xx — I yy )/I Z z the ellipticity, ijj the principal moment of inertia, 
B the projection of the magnetic dipole in the direction orthogonal to the rotation axis. The evolution 
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of the massive stars that give birth to pulsar being very fast, the rate can be derived directly from the 
star formation rate (see Eq. (11)). Considering the interval 8 — 40 M Q for the mass range of neutron 
star progenitors, and the initial mass function of Equation (12), Regimbau & Mandic (2008) found 
A - 10-2 Mg 1 . 

Normal radio pulsars, which are born with magnetic fields of the order of 10 12 — 10 13 G, rotational 
periods of the order of tens or hundreds of millisecond (Regimbau & de Freitas Pacheco 2000; Faucher- 
Giguere & Kaspi 2006; Soria et al. 2008), are not expected to contribute significantly to the GW signal 
(Regimbau & de Freitas Pacheco 2001b). However the population of newborn magnetars in which super- 
strong crustal magnetic fields (B ~ 10 14 — 10 16 G) may have been formed by dynamo action in a 
proto-neutron star with very small rotational period (of the order of 1 ms) (Duncan & Thompson 1992; 
Thompson & Duncan 1993), may produce a strong stochastic background in the frequency band of 
terrestrial detectors (Regimbau & de Freitas Pacheco 2006a). For these highly magnetized neutron stars, 
the distortion induced by the magnetic torque becomes significant, strongly enhancing the GW emission. 
In the case of a pure poloidal internal magnetic field matching to the dipolar field B in the exterior, the 
ellipticity is given by (Bonazzola & Gourgoulhon 1996; Konno et al. 2000): 

d8d2 

£ .=/W, (23) 



where j3 is a distortion parameter which depends on both the equation of state and the magnetic field ge- 
ometry. Using numerical simulations, Bonazzola & Gourgoulhon (1996) found that (3 can range between 
1 — 10 for a non-superconducting interior to 100 — 1000 for a type I superconductor and even take values 
larger than 1000 — 10 000 for a type II superconductor with counter rotating electric currents. Taking 
R = 10km for the radius, I zz — 10 45 gem 2 for the moment of inertia, and assuming that magnetars 
represent 10% of the population of NSs (Kouveliotou et al. 1998), we find that the stochastic signal is 
detectable with the Einstein Telescope after an observation time T = 1 yr and with a signal to noise 
ratio of 3 when | > 1.5 X 10~ 18 , giving the opportunity to put very interesting constraints on both B 
and j3. On the other hand, It has been suggested that the spindown could become purely gravitational if 
the internal magnetic field could is dominated by a very strong toroidal component (Cutler 2002; Stella 
et al. 2005), of the order of 10 16 G). In this saturation regime, the energy density increases as at low 
frequencies and reaches a maximum of il gw ~ 2.3 x 10 -8 around 760 Hz, giving a signal detectable by 
the Einstein Telescope with a signal-to-noise ratio of about 100 . 



5.3 Rotating Neutron Stars: Initial Instabilities 

5. 3. 1 Dynamical bar modes 

The gravitational stochastic background from core collapse supernovae could be enhanced by a number 
of proposed post-collapse emission mechanisms. One intriguing mechanism is the bar- mode dynam- 
ical instability associated with neutron star formation. These instabilities derive their name from the 
'bar-like' deformation they induce, transforming a disk-like body into an elongated bar that tumbles 
end-over-end. The highly non-axisymmetric structure resulting from a compact astrophysical object en- 
countering this instability makes such an object a potentially strong source of gravitational radiation and 
have been the subject of a number of numerical studies (Brown 2000; New et al. 2000; Shibata et al. 
2000; Saijo et al. 2001; Baiotti et al. 2007). Howell et al. have calculated the background resuting signal 
from this emission process using simulated energy spectra data, dE sw /du, from Shibata & Sekiguchi 
(2005), who performed the first three dimensional hydrodynamic simulations for stellar core collapse in 
full general relativity. Assuming a 20% occurrence of this instability, the authors find that the density 
parameter reaches a maximum of J7 gw ~ 4 x 10~ 10 around 600 Hz, and may be detectable with the 
Einstein Telescope with a signal to noise ratio of 3 after one year of integration. The optimistic event 
rate considered by Howell et al. is supported by suggestions that post collapse neutrino emission by the 
proto-neutron stars can induce contraction through cooling. This leads to increased spins though con- 
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servation of angular momentum (Shibata & Sekiguchi 2005). The implication here is that the instability 
can set in tens of milliseconds post collapse, increasing the rate of occurrence. 



5.3.2 r-modes 



The stochastic background from r-modes was first investigated by Owen et al. (1998) and then reviewed 
by Ferrari et al. (1999b). These estimates are based on the initial model of (Lindblom et al. 1998), 
which does not account for dissipation mechanisms such as the effect of the solid crust or the magnetic 
field, which may significantly reduce the gravitational instability. The spectral energy density of a single 
source is given by: 

^ = ^withz,G[0-Z4 up ], (24) 

dv I/2 up 

where i/ sup is 4/3 of the initial rotational frequency and Eq is the rotational energy lost within the 
instability window. For neutron stars with radius R = 10 km and mass M = 1.4 M Q the spectrum 
evolves as J7g W ~ 3 x 10 _12 £Vq where £ is the fraction of NS stars born near the keplerian velocity and 
which enter the instability window, until it reaches a maximum at 730,Hz. The Einstein Telescope may 
be able to detect this signal with a signal to noise ratio larger than > 3 for T = l yr if £ > 0.1%. One 
obtains similar constraints with the secular bar mode instability at the transition between Maclaureen 
and Dedekind configurations (Lai & Shapiro 1995). 



5.3.3 Collapse to quark matter 

It has been suggested that neutron stars could also undergo small core collapses after phase transitions, 
producing large amount of gravitational waves. Sigl (2006) calculated the background from phase transi- 
tion to quark matter in newly born NSs with millisecond periods, based on recent numerical simulations 
(Lin et al. 2006). Assuming that 1% of neutron stars are born fast enough to undergo the phase transi- 
tion, and that the energy released in the process represents about 5% of the rotational energy (~ 2 x 10 51 
erg), the energy density parameter may reach a maximum of fl gw ~ 10 -10 at kHz frequencies. 



5.4 Core Collapse Supernovas 

5.4.7 Core collapse supernovas to neutron stars 

After they have burnt all their nuclear combustible, massive stars may explode as type II supernovas. 
Their envelope is ejected while the core collapses to form a neutron star or a black hole (BH), de- 
pending on the initial mass of the progenitor, emitting a large amount of gravitational waves in the 
process. In a recent work, Howell et al. (2004) calculated the stochastic background that results from 
the birth of neutron stars at cosmological distances, using relativistic numerical models of core collapse 
(Dimmelmeier et al. 2002), and updating the previous study by Coward et al. (2001), based on newto- 
nian models (Zwerger & Mueller 1997). They considered three different GW waveforms, assumed to be 
representative of the three types of the catalog. Type I waveforms are characterized by a spike resulting 
from the core bounce followed by a ringdown, Type II by several distinct spikes and Type III shows 
large positive and smaller negative amplitudes just before and after bounce. In order to calculate the 
background spectrum, they assumed a flat cosmology, with O m = 0.3, — 0.7 and Hubble parameter 
Hq = 70 km s ~~ 1 Mpc - 1 (the 737 cosmology Rao et al. 2006), corresponding to the so-called concordant 
model derived from observations of distant type la supernovae (Perlmutter et al. 1999) and the power 
spectra of the cosmic microwave background fluctuations (Spergel et al. 2003), and considered three 
different models of the SFRs, finding no sensitive difference in the results. The NS progenitors were 
assumed to have masses between 8 — 25 M Q for a Salpeter IMF normalized between 0.1 — 125 M Q . The 
background is found continuous for Type II and rather a popcorn noise for Type I and III waveforms. 
The energy density parameter reaches a maximum of f2 gw ~ 3x 10~ 12 at around 700 Hz for Type I and 
n gw - 10~ 13 at 100 Hz and 800 Hz for Type II and Type III. 
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Besides the emission from the supernova bounce signal in the kHz range, it is expected that the 
large-scale convective overturn that develops in the delayed explosion scenario during the epoch of 
shock-wave stagnation, may emit a much stronger signal that may last for a few hundreds of ms before 
the actual explosion in the 1 Hz frequency range. Buonanno et al. (2005) estimated the background pro- 
duced by both ordinary supernovae and pop III stars using different numerical models of the GW wave- 
form (Fryer et al. 2004; Miiller et al. 2004). They showed that the signal is Gaussian below 1 Hz with an 
amplitude that may be at the level of the background expected from inflationary models. However, the 
authors stressed that these estimates remain uncertain by several orders of magnitude, essentially due to 
uncertainties on the parameters of the supernova GW emission. 



5.4.2 Core collapse supernovas to black hole 

The GW background from core collapse supernovas that result in the formation of black holes was first 
investigated by Ferrari et al. (1999a), using the relativistic numerical simulations of Stark & Piran (1985, 
1986) and later by de Araujo et al. (2002a) who found similar results assuming that all the energy goes 
into the ringdown of the I = m = 2 dominant quasi normal mode. For this mode, the frequency is given 
by (Echeverria 1989) 

v*(m,a) « " (25) 



with 



c 3 



A(a)= 0.63(1 -a) 03 ), (26) 

where M is the mass of the black hole, assumed to be a fraction a of the mass of the progenitor m, and 
a the dimensionless spin factor, ranging from for a Schwarzschild BH to 1 in the extreme Kerr limit. 
The spectral energy density has the simple expression: 

= eamc 2 8{u - vJM)) (27) 

av 

where e an efficiency coefficient. Previous numerical simulations of Stark & Piran (1985) gave an upper 
limit of e ~ 7 x 10~ 4 for an axisymmetric collapse, but accounting for more realistic scenarios, in 
particular the pressure reduction that triggers the collapse. Baiotti & Rezzolla (2006) obtained an effi- 
ciency of the order of 10~ 7 — 10~ 6 so 2 — 3 orders of magnitude smaller. Assuming that stars in the 
range 30 — 100 M can produce a BH, taking a — 10%, and a = 0.6, this simple model gives that the 
energy density ranges between 0.25 — 5.6kHz, with a maximum of Sl gw ~ e x 10~ 8 around 1650Hz, 
which means that an efficiency > 2 x 10~ 3 would give a signal detectable with a signal to noise ratio 
of 3 after one year of observation with the Einstein Telescope. Decreasing the minimal mass or a would 
narrow the spectrum and shift the maximum toward lower frequencies, while a change in the efficiency 
parameter e would only affect the amplitude. Increasing the spin factor or broadening its distribution 
broadens the spectrum and shifts the maximum toward larger frequency. Taking a = 20%, we find that 
the signal is detectable for efficiencies larger than 0.01%. 

In a recent work, Marassi et al. (2009) made used of the recent progress of numerical relativity, to 
review and extend the previous estimates of Ferrari et al. (1999a) for both population II and population 
III stars. The supernova rates were derived from the numerical simulations of Tornatore et al. (2007), 
which follows the star evolution, metal enrichment and energy deposition, and the GW signal from 
waveform derived from relativistic numerical simulations. The background is out of reach of the the 
first generation of detectors for Pop III stellar collapse, but could be detected by the Einstein Telescope 
for Pop II supernovas. Assuming 20 — 100 M Q for the mass range of BH progenitors, they found that the 
energy density reaches a maximum of fl gw ~4-7x 10~ 10 around 500 Hz for the model of Sekiguchi 
& Shibata (2005), giving a signal to noise ratio between 1.6 — 7.1 after one year of observation. In 
addition, they estimated the background from the collapse to neutron stars. Assuming 8 — 20 M Q for 
the mass of the progenitors, and the model of Ott (2005) which account for the g-mode excitation, they 
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Fig. 6 Energy density of the most promising astrophysical background contributions for 
ground based detectors, discussed in the text, magnetars (threshold detectable by ET), binary 
neutron stars, dynamical bar modes in proto neutron stars (courtesy of E. Howell) , r-modes 
assuming that 1% of newborn neutron stars cross the instability window, population II core 
collapse to neutrons (model of Ott et al. 2006) and to black holes (model D5a of Sekiguchi & 
Shibata 2005), courtesy of S.Marrasi) 



found that the energy density reaches a maximum f2 gw ~ 10~ 9 around 1000 Hz, giving a signal to noise 
ratio of 8.2. 

Similarly, Zhu, Howell & Blair (2010) estimated the GW signal created by all core collapse super- 
novae, to NS and BH, using Gaussian spectrum of the form 

^^ = Aexp{-(v-v*)/2o- 2 ) (28) 
dv 

shown to be a good approximations of the models of Ott et al. (2006). Based on simulated spectra 
of Dimmelmeier et al. (2008) and Sekiguchi & Shibata (2005), they considered different models with 
a ~ 500 and !/* = 200 — 800 Hz. They found that the signal may be detectable for efficiencies e > 10~ 5 
and e > 10~ 7 for Einstein Telescope. 

5.5 Capture by Supermassive Black Holes 

The emission from the various populations of compact binaries, which represent the main source of 
confusion noise for LISA, was studied intensively in the past decades (see for instance Kosenko & 
Postnov 1998; Ignatiev et al. 2001; Schneider et al. 2001; Farmer & Phinney 2002; Cooray 2004 for the 
extra-galactic contribution). The signal is expected to be largely dominated by white dwarf-white dwarf 
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(WD-WD), and in particular by the galactic population between 0.1 — 10 mHz (Yungelson et al. 2001; 
Nelemans et al. 2001; Benacquista et al. 2004; Edlund et al. 2005; Belczynski et al. 2005; Timpano 
et al. 2006; Benacquista & Holley-Bockelmann 2006), the extra-galactic contribution being one order 
of magnitude smaller (Kosenko & Postnov 1998; Schneider et al. 2001). In a recent paper, Barack & 
Cutler (2004) investigated the stochastic background created by unresolved captures by supermassive 
black holes (SMBHs) (Amaro-Seoane et al. 2007) . The capture rates for WDs, NSs and stellar BHs, 
which were extrapolated from the rates derived by Freitag (2003) for our galaxy, represent the main 
source of uncertainties, ranging between 4 x 10~ 8 — 4 x 10 -6 M^ 8 yr _1 for WD captures and between 
6 x 10~ 8 - 6 x 10~ 7 Afg /8 yr^ 1 for NS and BH captures, M 6 being the mass of the SMBH in units of 
10 6 M©. In Figure , the most optimistic and pessimistic models are compared to the LISA instrumental 
noise and to the WD-WD galactic foreground derived by Bender & Hils (1997). For the most optimistic 
rates, the resulting background may contribute to the LISA confusion noise, raising the effective LISA'S 
overall noise level by a factor of ~ 2 in the range 1 — 10 mHz, where LISA is most sensitive. 




1E-4 1E-3 0.01 



f (Hz) 



Fig. 7 gravitational strain in Hz -1 / 2 , corresponding to optimistic (grey continuous curve) and 
pessimistic (grey dashed curve) compact object captures (Barack & Cutler 2004), along with 
the LISA instrumental noise (black) and the WD-WD foreground (black). 



6 CONCLUSIONS AND FUTURE PROSPECTS 

Gravitational Waves ground based experiments, after a decade of detector installation and commission- 
ing, have reached or surpassed their design sensitivities, opening a new window into the Universe, The 
first generation LIGO interferometers have already put interesting astrophysical constraints on the ellip- 
ticity of the Crab pulsar (below the spindown limit). With advanced detectors, we expect to see at least 
close compact binary coalescences, while the third generation detector Einstein Telescope and the space 
detector LISA should bring GW astronomy to the next level, when it is possible to address a range of 
problems on a wide variety of astrophysical sources but also in fundamental physics and cosmology. 
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The cosmological stochastic background is often seen as the Graal of GW astronomy since it would 
give a snapshot of the very early stages of the Universe, up to a fraction of second after the Big Bang. The 
astrophysical background is also promising since it would provide information on the physical properties 
of compact objects and their evolution with redshift, such as the mass of neutron stars or black holes, 
the ellipticity and the magnetic field of neutron stars, the angular momentum of black holes, the rate of 
compact binaries. We have shown in the previous sections that astrophysical models are out of reach 
of the first generation of detectors but with advanced detectors, and particularly the third generation 
Einstein Telescope, upper limits could put very interesting constraints on the equation of state and the 
magnetic field of magnetars, the distribution of the birth rotation period of newborn neutron stars and 
models of core collapse supernovas. Moreover, unless we overestimate the rate by orders of magnitude, 
we should be able to see the background from coalescing double neutron star binaries. 

On the other hand, the astrophysical contribution may be a noise masking the cosmological back- 
ground, and also a confusion foreground where the detection of individual high redshift standard candles 
needed to infer dark energy may become difficult (Regimbau & Hughes 2009). In this context, modeling 
the astrophysical background as precisely as possible to extract informations on its strength, frequency 
range and statistical properties, anything that may help distinguish it from the cosmological signal or 
separate overlapping sources is crucial. 

Another important task in the next few years will be to adapt our actual methods to the triangle 
configurations of ET and the spatial interferometer LISA. It is not sure yet that we can ever get rid of 
the correlated noise, but we can certainly be able to reduce it by correlating interferometers formed by 
specific combinations of the three arms, with the extra complication for LISA that the triangle is moving. 

Acknowledgements The author thanks Vuk Mandic for providing the landscape plot for cosmological 
background, Eric Howell and Stefania Marassi for their contribution to the bar mode instability and the 
core collapse supernovas sections. 
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